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Abstract 

Although neutrons do not couple directly to the superconducting order parameter, they have nevertheless 
played an important role in advancing our understanding of the pairing mechanism and the symmetry of 
the superconducting energy gap in the iron arsenide compounds. Measurements of the spin and lattice 
dynamics have been performed on non-superconducting 'parent' compounds based on the LaFeAsO ('1111') 
and BaFe2As2 ('122') crystal structures, and on electron and hole-doped superconducting compounds, us- 
ing both polycrystalline and single crystal samples. Neutron measurements of the phonon density-of-state, 
subsequently supported by single crystal inelastic x-ray scattering, are in good agreement with ab initio 
calculations, provided the magnetism of the iron atoms is taken into account. However, when combined 
with estimates of the electron-phonon coupling, the predicted superconducting transition temperatures are 
less than 1 K, making a conventional phononic mechanism for superconductivity highly unlikely. Measure- 
ments of the spin dynamics within the spin density wave phase of the parent compounds show evidence of 
strongly dispersive spin waves with exchange interactions consistent with the observed magnetic order and 
a large anisotropy gap. Antiferromagnetic fluctuations persist in the normal phase of the superconducting 
compounds, but they are more diffuse. Below Tc, there is evidence in three '122' compounds that these 
fluctuations condense into a resonant spin excitation at the antiferromagnetic wavevector with an energy 
that scales with Tc. Such resonances have been observed in the high-Tc copper oxides and a number of 
heavy fermion superconductors, where they are considered to be evidence of d-wave symmetry. In the iron 
arsenides, they also provide evidence of unconventional superconductivity, but a comparison with ARPES 
and other measurements, which indicate that the gaps are isotropic, suggests that the symmetry is more 
likely to be extended-s-t wave in character. 
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1. Introduction 

Since the discovery of superconductivity in iron 
arsenide compounds [H, 13,13' neutron scattering 
experiments have made significant contributions to 
our understanding of the underlying physics. Early 
neutron diffraction results generated considerable 
excitement because they revealed remarkable simi- 
larities with the high-temperature copper oxide su- 
perconductors. For example, in both the iron ar- 
senides and the cuprates, superconductivity arises 
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when an antiferromagnetically ordered phase has 
been suppressed by chemical doping [^. Neutron 
scattering continues to be essential in determin- 
ing the magnetic and structural phase diagrams of 
these materials as a function of dopant concentra- 
tion or applied pressure On the other hand, 
neutrons have also identified important differences 
with the cuprates, such as the reduced size of the 
ordered moments and their extreme sensitivity to 
structural modifications [1, [E I3l- Elastic neutron 
scattering, which probes static magnetic and struc- 
tural correlations, is discussed in more detail in an- 
other article in this issue Q. The purpose of this 
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review is to summarize the results of inelastic neu- 
tron scattering, which probes dynamic correlations 
involving phonons and spin fluctuations, both of 
which are candidates for binding the superconduct- 
ing electron pairs. 

With the discovery of any new family of super- 
conductors, the first task is to determine whether 
the critical temperature can be explained by 
electron-phonon coupling within a conventional 
BCS theory. This question is usually addressed 
within the formalism of Eliashberg theory 8], in 
which the superconducting energy gap is expressed 
in terms of a spectral density function derived from 
the phonon density-of-states (PDOS) weighted by 
electron-phonon matrix elements. We will review 
inelastic neutron scattering, mostly on polycrys- 
talline samples, that have been used to estimate the 
PDOS [9'] and validate the results of first principles 
density functional calculations 0, Ei. In broad 
terms, the agreement between theory and experi- 
ment is very good, although some modes are ex- 
tremely sensitive to the spin state assumed in the 
theoretical estimates. Although there are subtle, 
so far unexplained, anomalies that will require fur- 
ther single crystal measurements to resolve fully, 
the early consensus is that the electron-phonon cou- 
pling is too weak by a factor of about five to explain 
the observed critical temperatures. 

If phonons are not responsible for the super- 
conductivity, spin fluctuations offer an alternative 
bosonic spectrum to mediate the electron pair- 
ing. As flrst shown by de la Cruz et al 4], the 
ground state of the non-superconducting parent 
compounds is a spin density wave, whose transition 
occurs close to a tetragonal-orthorhombic struc- 
tural transition. Both chemical doping and, in con- 
trast to the cuprates, pressure [I4I can be used to 
suppress antiferromagnetic order and induce super- 
conductivity. One of the key questions to resolve is 
what drives the magnetism. Band structure cal- 
culations show that hole pockets at the F-point 
and the electron pockets at the M-point can show 
strong nesting with a sharp peak in the Lindhard 
susceptibility at Q=(7r,7r), using tetragonal nota- 
tion [13|. Since this is also the wavevector of mag- 
netic order, it is natural to propose that the un- 
doped arsenides are itinerant spin density waves 
like chromium, which would provide an explana- 
tion for the reduced size of the ordered moment. 
However, some have argued that the pnictides are 
in fact close to a Mott insulating phase and that 
the reduced moments are due to frustration caused 



by competing superexchange interactions IJ, |15| , 
so the strength of electron correlations remains an 
important issue. Inelastic neutron scattering exper- 
iments have shown that the spin waves are strongly 
dispersive, with velocities that are not inconsistent 
with an itinerant spin density wave, and more three- 
dimensional than the cuprates. They also show sub- 
stantial energy gaps that are not fully understood. 

Finally, we review three reports of a spin reso- 
nant excitation seen so far only in the '122' com 
pounds [F" 
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With chemical doping, the 
spin fluctuations become more diffuse in the nor- 
mal state, though still centered at the antiferro- 
magnetic wavevectors. However, below Tc, these 
fluctuations condense into a resonant excitation 
that is localized in both momentum transfer, Q, 
and energy transfer, uj. Such resonant excita- 
tions have been observed in a wide range of high- 
temperature copper oxide superconductors [19] as 
well as, more rece ntly , several heavy fermion su- 
perconductors El IS El , where they are con- 
sidered to be evidence of d-wave superconductivity 
2J,|25|. In the case of the iron arsenides, a compar- 
ison with ARPES data suggests that the resonance 
is evidence of extended-s± wave symmetry [2H [27l|. 
in which the disconnected hole and electron pockets 
have energy gaps of opposite sign. Since there is no 
angular anisotropy of the gap in this symmetry, it is 
difficult to verify by the techniques used in the cop- 
per oxides. Inelastic neutron scattering is so far the 
only probe that has provided phase-sensitive infor- 
mation about the unconventional symmetry of the 
energy gap in the iron arsenide superconductors. 

2. Lattice Dynamics 

Within Eliashberg theory, the superconducting 
gap equation is expressed in terms of an electron- 
jhonon coupling, A, and a Coulomb repulsion, /i* 
^. The electron-phonon coupling is derived from 
the spectral density function a^F{u}), in which the 
phonon modes, represented by the bare phonon 
density-of-states (PDOS), F{uu), are weighted by 
electron-phonon matrix elements. In many super- 
conductors, there is little difference in the func- 
tional forms of q;^F(w) and F{uj), i.e., the phonons 
are all coupled to the electronic states equally 
strongly, but there are exceptions. For example, 
in MgB2, the electron-phonon coupling of the E2g 
modes is particularly strong because their energies 
are governed by strongly covalent boron-boron a- 
bonds. In principle, a'^F{uj) can be determined 



by inverting tunneling data, but this is not al- 
ways available so neutron scattering measurements 
of i^(a;) play a valuable role in providing initial esti- 
mates of Tc and in validating ab initio calculations. 

In BCS superconductors, numerical calculations 
based on the Eliashberg equations have shown that 
the critical temperature can often be approximated 
using the Allen-Dynes equation [2^ 
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where uin is a logarithmic phonon average defined 
in equation 2.16 of ref. ,8], ^* is the Coulomb repul- 
sion, and the electron-phonon coupling, A, is given 

by 
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The Coulomb repulsion, 11* , is normally treated as 
a phenomenological parameter but it is possible to 
make accurate predictions of the critical tempera- 
ture from first principles calculations without any 
adjustable parameters (29| . 

There were two ab initio calculations of the 
phonon density-of-states in LaFeAsO before the 
first neutron measurements were published, and 
both are in broad agreement with each other 11, 
[loj . Singh and Du estabhshed the main features 
of the band structure in the Local Spin Density 
Approximation, showing that the low-lying elec- 
tronic states arise from five iron d-bands spanning 
an energy range of -2.1 eV to 2.0 eV, with the oxy- 
gen and arsenic p-states well below the Fermi level 
[lO| . The Fermi surface consists of two-dimensional 
cylinders, with two hole pockets at the F-point and 
two electron pockets at the M-point, derived mainly 
from dxz and dyz states, along with a fifth more 
three-dimensional hole pocket, centered at Z, de- 
rived from dj,2 states hybridized with As p-states 
and La orbitals. 

The resulting phonon excitations are spread over 
70 meV, with those over 40 meV representing oxy- 
gen modes. Those at lower frequency are of mixed 
lanthanum, iron, and arsenic character. As sev- 
eral of the electron bands are two-dimensional, it 
is not surprising that the optic modes show little 
dispersion along the F-Z direction, but the acoustic 
phonons are fairly three-dimensional in character 
with a Debye temperature of 340 K. The calculated 
phonon density-of-states has three main peaks at 
approximately 12, 21, and 34meV. 




Figure 1: Inelastic neutron scattering from 

LaFeAsO0.89F0.11 as a function of momentum trans- 
fer, Q, and energy transfer, AE (denoted as uj in the 
text), measured at 10 K with an incident neutron energy of 
100 meV . The units of the intensity scale are arbitrary. 



On the basis of such band structure calculations, 
Boeri et al calculate the electron-phonon coupling, 
which they find to be evenly distributed among all 
the modes [ll|. When they integrate this coupling 
over all frequencies, they derive A — 0.21, which 



is much smaller than any other known electron- 
phonon superconductors. In the Allen-Dynes ap- 
proximation, they estimate Tc = 0.5 K, from their 
calculated value of win = 205 K, assuming fi* = 0. 
A more accurate numerical Migdal-Eliashberg cal- 
culation only increases this to 0.8 K, and they esti- 
mate that A would have to be a factor five stronger 
to generate the observed Tc of 26 K. Their calcula- 
tions are based on the non-superconducting parent 
compound, LaFeAsO, but adding electrons through 
fluorine doping would, in the rigid band approx- 
imation, reduce the electronic density-of-states at 
the Fermi level, and so tend to reduce Tc even fur- 
ther. Orbitals that might be expected to produce 
an enhanced electron-phonon coupling, such as the 
dr^i_yi orbitals that are directed along the Fe-Fe 
bond, are too far from the Fermi level to have a 
strong influence. 

Even before any measurements were reported, it 
therefore seemed unlikely that a conventional BCS 
mechanism could explain the elevated transition 
temperatures in the iron arsenides. Nevertheless, 
it was important to validate these predictions from 
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Figure 2: Generalized phonon density of states, G{lj), of 
LaFeAsOi— a;F2: measured with an incident neutron energy 
of 130 meV and normalized to an integral of 1.0 after cor- 
rection for the Bose population and Debye- Waller factors, 
detector efficiency, and multiphonon scattering Q. (a) Ex- 
perimental G{u!) measured at 35 and 300 K for x = 0.1. 
(b) Experimental G{uj) for x = and x = O.l measured at 
300 K. (c) First-principles calculation of the phonon density- 
of-states of LaFeAsO based on the band structure of Singh 
and Du O. 

experiment. Although inelastic x-ray scattering can 
also be used to measure phonon dispersion relations 
(and has been used in the iron pnictides 13, 3l|). 
neutron scattering is an extremely efficient method 
of determining the phonon density-of-states. In a 
multicomponent system, the inelastic neutron scat- 
tering law in polycrystalline samples is given by 
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where ct^ and Mi are the neutron scattering cross 
section and atomic mass of the z*'' atom and n{u}) = 
[exp {huj/kBT) — l]~^ istheBosepopulation factor . 
The generalized PDOS, G(w) = E»G,(a;), where 



Gi{uj) is defined as 



which, in turn, defines the Debye-Waller factor, 
W,(Q) through 



2M, 



duj- 



-[2n{uj) + l] (5) 



a;(j, q) and ei(j, q) are the frequencies and eigen- 
vectors, respectively, of the phonon modes. 

Inelastic neutron scattering therefore measures 
a sum of the partial phonon density-of-states of 
each constituent element weighted by ai/Mi. In 
LaFeAsO, these weights, relative to the value for 
oxygen, are 0.23, 0.76, and 0.27, for La, Fe, and 
As, respectively [9|, so the low energy scattering is 
dominated by iron modes. The generalized PDOS 
also differs from the base PDOS, F(a;), because of 
additional weighting by the eigenvectors (see Equa- 
tion[3]) . Osborn et al argued that this makes the use 
of G(w) preferable to F(a;) as an approximation to 
a^F(a;), because these eigenvectors also enter into 



electron-phonon matrix elements 33] ■ Neverthe- 
less, the safest way to compare ah initio calculations 
to the neutron data is to calculate the generalized 
PDOS directly as was done by Bohnen et al in the 
case of MgBa [H. 

Strictly speaking. Equation [3] only applies to ma- 
terials in which the neutron cross sections of the 
constituent elements, di, are entirely incoherent, 
whereas both iron and arsenic have strongly coher- 
ent cross sections. This means that there will be 
strong deviations from the simple dependence 
of the scattering intensity caused by variations in 
the eigenvectors within each Brillouin zone. Nev- 
ertheless, if the data are taken on polycrystalline 
samples, so that the measured scattering is spher- 
ically averaged, and integrated over a sufficiently 
broad range of Q, they can be used to derive a good 
approximation to the generalized PDOS even when 
the tTi are coherent. This is known as the incoher- 
ent approximation, which must be satisfied for the 
analysis to be reliable. In MgB2, for example, there 
were reports of a low-energy peak that were then 
shown to be an artifact of insufficient Q-averaging 
[s^ ]. However, most of the reported PDOS mea- 
surements in the iron arsenides were taken on neu- 
tron spectrometers using relatively high incident 
energies and summed over a large range of momen- 
tum transfers, ensuring that the conditions for the 
incoherent approximation are met. 




Figure 3: The experimental phonon spectra of Sro.6Ko.4Fe2As2 and Cao.6Nao.4Fe2As2, measured in neutron energy gain with 
an incident neutron energy of 3.1 meV ^5] i compared to BaFe2 As2, measured in energy loss with an incident energy of 57.5 meV 
[3^ . All the phonon spectra are normalized to unity. 



The first reported phonon measurements in the 
'1111' system covered a Umited energy range up 
to 20meV [sl], but subsequent experiments have 
covered the entire phonon spectrum . Christian- 
son et al used samples of both non-superconducting 
LaFeAsO and superconducting LaFeAsOi_a;Fa;, 
with X ~ 0.1 prepared by two different synthesis 
groups, based at Oak Ridge National Laboratory 
and Ames Laboratory. The data were collected on 
the newly commissioned Fermi chopper spectrom- 
eter, ARCS, at the Spallation Neutron Source in 
Oak Ridge. Using incident neutron energies of 130, 
60, and 30meV, supplemented by some triple-axis 
measurements, they observed peaks at 12, 25, 31, 
40, and 60meV. These features seem to be unaf- 
fected either by temperature reduction, apart from 
a sharpening of the peaks, or by dopant concen- 
tration. In particular, there does not appear to be 
any strong phonon renormalization either at the or- 
thorhombic transition in the non-superconducting 
samples, or at Tc in the superconducting samples. 

A comparison of their data with the ab initio cal- 
culations showed good qualitative agreement in the 
overall energy scale of the phonons and the energies 
of most of the peaks. The most notable discrep- 
ancy was in the location of the 31 meV peak, which 
is distinctly softer than the theoretical prediction. 
A similar discrepancy was noted in inelastic x-ray 
scattering, where it was attributed it to a 30% re- 



duction in the Fe-As force constants 30] . The length 
of the Fe-As bond shows a remarkable sensitivity 
to the iron spin state; calculations that underesti- 
mate the iron magnetism also substantially under- 
estimate the bond lengths, and therefore overesti- 
mate the bond energies. The fact that there is such 
little variation in the phonon peak energies is evi- 
dence therefore that, locally, the iron spin state is 
remarkably robust, persisting above the SDW tran- 
sition temperature in LaFeAsO and surviving the 
destruction of SDW order in LaFeAsOo.gFo.i. This 
is borne out by measurements of the spin dynamics 
discussed in the next section. 

There have also been phonon measurements 
performed on the '122' compounds, firstly with 
experiments on non-superconducting BaFe2As2 
[34| and then on superconducting Sro.6Ko.4Fe2As2 
(Tc=32K) and Cao.6Nao.4Fe2As2 (Tc=21K) 
The measurements were taken on the IN4 spec- 
trometer, using an incident neutron energy of about 
60meV, and the IN6 spectrometer, using an inci- 
dent energy of 3.1 meV, both at the Institut Laue 
Langevin, France. The IN6 data were measured 
in neutron energy gain, which requires an elevated 
temperature. With the absence of light oxygen 
ions in this structure, the maximum phonon energy 
is just under 40meV. However, in other respects, 
the measured phonon spectra of the different '122' 
compounds look similar to the '1111' compounds 



and to each other. There are peaks at 12, 22, 27, 
and 34meV in the barium and strontium samples 
but, in the calcium sample, the 22meV peak ap- 
pears to have shifted down to below 18meV. Since 
the bond lengths are shorter in Cao.6Nao.4Fe2As2, 
which would normally increase the mode energies, 
there must be some change in the bonding char- 
acteristics. Mittal et al note that there is a slight 
stiffening from 300 K to 140 K of the higher energy 
peaks in the IN6 data, but a softening of the acous- 
tic modes, which they suggest is a sign of electron- 
phonon coupling [sB]. 

In summary, neutron scattering studies of the lat- 
tice dynamics of the iron arsenides are broadly con- 
sistent with ah initio calculations, although there 
are some unexplained anomalies that will require 
more detailed single crystal measurements before 
they are explained satisfactorily. The integrated 
electron-phonon coupling is estimated to be far too 
weak to be responsible for the superconductivity. 
However, the sensitivity of the Fe-As bond length 
to the iron spin state shows that there are potential 
sources of electron-phonon coupling that may need 
to be investigated more fully before some contribu- 
tion from electron-phonon coupling is definitively 
ruled out. 



3. Spin Dynamics 

The non-superconducting parent compounds, 
such as LaFeAsO or BaFe2As2, undergo twophase 
transition with decreasing temperature 0, 0, [37j . 
The first is a structural transition from the high- 
temperature tetragonal phase to a low-temperature 
orthorhombic phase, which is closely followed by (or 
sometimes coincident with) a second magnetic tran- 
sition [5|. The low-temperature antiferromagnetic 
structure is not the checkerboard order that would 
be expected from nearest-neighbor antiferromag- 
netic interactions, but a stripe phase, which implies 
the presence of competing interactions (see Fig. 4). 
Yildirim has shown that this stripe structure re- 
sults from an inherent frustration produced by a 
strong next-nearest-neighbor antiferromagnetic ex- 
change between iron spins on the square planar 
lattice [il]. The observed structure is stable pro- 
vided the nearest-neighbor exchange, Ji, along the 
square edges, and next-nearest-neighbor exchange, 
J2, along the square diagonals, satisfy, J2 > Ji/2. 
He postulates that the structural transition is a 
means of relieving this frustration, which implies 
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Figure 4: The crystal and magnetic structure of BaFe2As2 
(taken from Ref. The unit cell contains two layers 

of Fe2As2 tetrahedra (Fe, blue spheres; As, yellow spheres), 
separated by planes of barium atoms (red spheres) . The blue 
arrows show the observed ordering of the iron spins. The red 
arrow shows the spacing of the antiferromagnetic stripes. 



the existence of strong short-range spin correla- 
tions above the antiferromagnetic transition tem- 
perature, TsDWi consistent with the discussion in 
the previous section. 

As stated in the introduction, there are two al- 
ternative explanations for the origin of the mag- 
netic interactions. The first is that the antifer- 
romagnetism is produced by a nesting instability 
coupling the F-centered hole pockets and the M- 
centered electron pockets. This assumes an itin- 
erant picture of weakly interacting electrons, that 
is consistent with the reduced size of the mag- 
netic moments (0.36/iB in LaFeAsO 4] or 0.87/^s in 
BaFe2As2 [13]), which can be explained by density 
functional theory provided the experimental lattice 
parameters are used [39] . However, there are alter- 
native models, which assume much stronger elec- 
tron correlations, with moments reduced by frus- 
tration [14, ,15.]. 

Inelastic neutron scattering is the most direct 
way of determining the spin wave excitations of the 
antiferromagnetically ordered compounds. What- 
ever the origin of the magnetic interactions, it is 
usually possible to analyze the measured disper- 
sion using a Heisenberg model including single-ion 
anisotropy terms that are required to produce the 
observed energy gaps 4^, 41] 



H — ^ Si • Sj + 



(6) 



Ewings et al derive explicit solutions of this 
Hamiltonian, giving both the energies and spin 
wave cross sections [i^, but, at low energies, the 
spin wave dispersion can be approximated by [i^ 



(7) 



where q is the reduced wavevector relative to the 
antiferromagnetic zone center, A is the anisotropy 
gap, and Vxy and Vz are the in-plane and c-axis spin 
wave velocities. 

Spin wave measurements have been reported in 
a number of the '122' compounds, SrFe2As2 [i^ ]. 
CaFezAsa [11,113, and BaFesAss 0,111]. Most of 
these are single crystal measurements using triple- 
axis spectrometers, although Ewings et al and 
Diallo et al [i^l used pulsed neutron source Fermi 
chopper spectrometers, MERLIN and MAPS re- 
spectively, at ISIS, to extend the energy range to 
100 meV and above. Broadly speaking, all these ex- 
periments are consistent. They observe extremely 
steep spin waves emerging from the SDW wavevec- 
tor with a substantial energy gap, ranging from 
6.9 meV in SrFeaAsa Q to 9.8 meV in BaFeaAsa 
[45| . The value of the spin wave velocity is more 
uncertain because it is difficult to resolve the prop- 
agating modes in constant energy scans when the 
dispersions are so steep, except at very high en- 



ergy [4j|. However, estimates can be derived by 
performing model simulations including the instru- 
mental resolution. The in-plane estimates vary 
from 280mcVA in BaFcaAsa 0, Si] to 420meVA 
in CaFeaAsa (42|. Low-energy measurements pre- 
dicted zone boundary energies of about 175 meV 
(iol . list , which is consistent with the first single 



crystal measurements using pulsed neutrons [4J|, 
although the spin waves above 100 meV appear to 
be heavily damped. The out-of-plane velocities are 
smaller but still substantial, reduced from the in- 
plane values by a factor of 2 42 1 to 5 [il], so 



the magnetic order is truly three-dimensional. The 
temperature dependence of the spin wave scattering 
is consistent with Bose statistics, with the energy 
gap renormalizing to zero at T^sdw- However, there 



120 



Log(lntensity) 




25 



20 



Intensity 
(mbarn/sr/meV/f.u.) 



9 

CD 

E 15 



LU 



10 



1 

-(b) 


it 











0.5 



Q (A-^) 



1.5 



n 



Figure 5: Neutron scattering spectra from a polycrystallinc 
sample of BaFe2As2 at 7K [4011 . The data were taken on the 
MERLIN (ISIS) spectrometer using incident neutron ener- 
gies of (a) 200 meV and (b) 50 meV. The pillars of scatter- 
ing show the steep spin wave excitations emerging from the 
and (1,1, positions (using the tetragonal Brillouin 
zone) at Q=1.2A and 2.6 A, respectively. 



is evidence of the persistence of two-dimensional 
short-range spin correlations in BaFeaAsa above 
TsDWi in the form of q uasielastic scattering cen- 
tered on (^,^,0 rods |45| . 

If the Heisenberg model includes next-nearest- 
neighbors, as required to stabilize the stripe struc- 
ture, there are four exchange constants: J\a and J\h 
couple antiferromagnetic and ferromagnetic near- 
est neighbors, respectively; J2 couples next-nearest 
neighbors within the plane; and Jz is the out- 
of-planc coupling. In SrFeaAsa, Zhao et al esti- 
mated the exchange interactions to be J\a + 2 J2 ^ 
100 meV, with ~ 5meV (i^. These values 
are consistent with McQueeney et aVs analysis of 




an unusual temperature dependence of the excita- 
tions in LaFeAsO is also consistent with itinerant 
spin density waves [48| . 

In summary, the neutron measurements in the 
parent compounds of the iron arsenides provide a 
consistent picture of spin wave dispersions, that are 
extremely steep and anisotropic, although clearly 
three-dimensional, with a very large energy gap. 
The exchange parameters are consistent with the 
observed stripe phase, and are not inconsistent with 
itinerant SDW models, although it is too early to 
rule out more localized magnetic models based on 
strong but frustrated superexchange interactions. 



Figure 6: Constant energy scans performed on triple-axis 
spectrometers on SrFe2As2 at 160 K, showing the broadening 
as a function of increasing energy resulting from the spin 
wave dispersion [4^ . 



CaFe2As2, where they estimate Jia = 41meV, Ju 
= lOmeV, J2 = 21meV, and = 3meV. Note 
that the value of J2 is sufficiently large to stabilize 
stripe antiferromagnetism. 

In an itinerant SDW, the spin wave velocity pro- 
duced by the excitation of particle-hole pairs across 
the nested surfaces is given by y^v^vt, where Ve 
and Vh are the electron and hole velocities of their 
respective bands [i^. In LaFeAsO, band struc- 
ture calculations estimate that the in-plane elec- 
tron velocity is a factor of 3 greater than the in- 
plane hole velocity, which is, in turn, a factor of 2.5 
greater than both the c-axis hole and electron veloc- 
ities The absolute values are uncertain since 
the calculated bandwidths are greater than ARPES 
measurements [13]. They do, however, give a ratio 
of in-plane to c-axis spin wave velocities of about 
4, which is within the experimental range. Zhao et 
al go further and show that the absolute values are 
in reasonable agreement with band structure, after 
renormalizing the bandwidths [43| . 

The damping of the high-energy excitations seen 
in CaFe2As2 could be due to the decay of spin waves 
into electron-hole pairs within a Stoner continuum 
(4^ . More speculatively, Matan et al report the 
existence of excess scattering above 20 meV, which 
is more two-dimensional than the spin waves, and 
which they suggest is consistent with the existence 
of excess scattering in the Stoner continuum. How- 
ever, the data from Ewings et al do not show evi- 
dence for this anomaly 40] , so it is difficult to draw 
any firm conclusions. Finally, there are claims that 



4. Resonant Spin Excitations 

The previous two sections concerned inelastic 
neutron scattering results that provide the context 
within which superconductivity develops in the iron 
arsenides, but do not directly address the super- 
conducting state itself. That has been the tra- 
ditional role for neutron scattering since, in con- 
ventional superconductors, neutrons do not couple 
directly to the superconducting order parameter. 
However, neutron scattering results on other un- 
conventional superconductors over the past twenty 
years has shown that they can provide direct in- 
formation concerning the superconducting energy 
gap, in particular, its symmetry. Surprisingly, in 
the case of the iron arsenides, it is currently the 
only technique to give phase-sensitive evidence of 
the gap symmetry, since many of the techniques 
commonly used in other superconductors do not ap- 
ply- 

Inelastic neutron scattering measures directly the 
dynamic magnetic susceptibility of the band elec- 
trons. Normally, this signal is too diffuse for useful 
measurements, but there are cases where the mag- 
netic response is strongly enhanced at particular 
energies and wavevectors making such experiments 
feasible. The most common example is in transi- 
tion metal magnets such as iron and cobalt, where 
the RPA susceptibility of the itinerant d-electrons 
has poles corresponding to propagating spin wave 
modes [i^. In an itinerant SDW model, the spin 
waves described in the previous section fall into this 
category. No such divergences exist in the dynamic 
magnetic susceptibility of conventional s-wave su- 
perconductors, but it turns out that they can exist 
when the superconducting energy gap changes sign, 
either within a single Fermi surface or between two 
disconnected Fermi surfaces. 




Figure 7: Inelastic neutron scattering on Bao.6Ko.4Fe2As2, measured using an incident neutron energy of 60meV below and 
above Tc, i.e., 7K (a) and 50K (b) respectively, and 15meV below and above Tc, (c) and (d) respectively. The data shows 
the transfer of spectral weight from diffuse spin fluctuations centered at 1.15 in the normal state into a resonant spin 
excitation at an energy transfer of 15meV in the superconducting state. The color scale is in units of mbarns/sr/meV/mol 




Specifically, there is a strong enhancement of the 
dynamic magnetic susceptibility due to coherence 
factors that appear because of the anomalous Green 
function of the superconducting phase. The follow- 
ing term appears in the non-interacting susceptibil- 
ity [50,[5li[2i,|23. 



1 



^k+qCk + AkAk+q 



(8) 



where the energy of the superconducting quasipar- 
ticles, £^k 



-r A^. Here, ^k are the single 
electron energies of the normal-state and Ak are 
the values of the energy gap at points k on the 
Fermi surface. 

Equation [S] defines the principal characteristics 
of the resonant spin excitations. If Ak+Q — Ak, 



as in a conventional s-wave superconductor, Equa- 
tion 8 vanishes on the Fermi surface (^k = 0)- In 
this case, there is no pole when this non-interacting 
susceptibility is introduced into an RPA expres- 
sion for the interacting susceptibility. However, if 
Ak+Q = — Ak, i.e. when Q connects two Fermi 
surface points whose superconducting order param- 
eters have opposite sign, then Equation 8 is maxi- 
mal, resulting in a pole at Q for some energy less 
than 2 A. 

Such resonant excitations are now believed to be 
a universal feature of the high-temperature copper 
oxide superconductors with observations in a 
large number of different systems [H, [HI, [5l|, [HI, 
55l . [56| . The resonance energy scales approximately 
as (jjQ ~ SfcsTc. They are commonly taken as evi- 
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Figure 8: The inelastic neutron scattering from 
Bao.6Ko.4Fe2As2 integrated over a Q-range of 1.0 to 
1.3 and an range of 12.5 to 17.5 meV. The integration 
range corresponds to the region of ma^cimum intensity of 
the resonant excitation observed below Tc (see Fig. 7). The 
dashed line is a guide to the eye below Tc and shows the 
average value of the integrals above Tc [tTj . 



dence of d.j.2_y2 symmetry, in which the energy gap 
changes sign within a single Fermi surface. Strik- 
ingly, there are now several reports that similar res- 
onant excitations are present in heavy fermion su- 
perconductors 2^ 21, 22, 2^], where Tc can be as 
low as 0.7 K. Nevertheless, the energy of the reso- 
nance appears to obey the same scaling relations, 
with cl;o/2A, where A is the maximum value of the 
superconducting energy gap, taking values between 
0.62 and 0.74, a remarkable observation given that 
Tc varies by over two orders of magnitude. 

Figure 7 shows inelastic neutron scattering data 
taken on the MERLIN spectrometer, at the ISIS 
Pulsed Neutron Source. The polycrystalline sam- 
ple of Bao.6Ko.4Fe2As2 used in these measurements 
was optimally doped with a Tc of 38 K. The intense 
scattering at high Q and low energy results from 
phonon and elastic nuclear scattering, respectively, 
but the remaining scattering results from spin fluc- 
tuations. In the normal phase, there is an echo of 
the sharp spin waves seen in Fig. 5, with a pillar 
of scattering at the antiferromagnetic wavevector, 
Q=1.15A . However, the scattering is considerably 
more diffuse representing short-range spin fluctua- 
tions rather than propagating spin waves. It also 
persists down to lower energy. As the temperature 
is lowered through Tc, there is a transfer of spec- 
tral weight into the resonant spin excitation at an 
energy transfer of 15 meV. Since the maximum gap 
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Figure 9: (a) Calculated imaginary part of the RPA spin 
susceptibility at the SDW wave vector Qafm a-s a function 
of frequency in the normal and superconducting states. The 
red, dotted blue, and solid blue curves correspond to the 
total RPA susceptibility for non-superconducting, d^2_y2 
symmetry and extended-s± symmetry models.. The thin 
(black) curves refer to the partial RPA contributions for the 
interband and intraband transitions in the s± superconduct- 
ing state, (b) Calculated imaginary part of the total RPA 
spin susceptibility in the s± state as a function of frequency 
and momentum along Q={h,h) [ssll . 



seen in ARPES data is 12meV [57|, a;o/2A - 0.58 
in good agreement with the scaling in other uncon- 
ventional superconductors. The temperature de- 
pendence of this resonant spin excitation is shown 
in Fig. 8, where it behaves like an order param- 
eter, as also observed in the other unconventional 
superconductors. 

The observation of a spin resonance does not nec- 
essarily imply d-wave symmetry. In fact, ARPES 
data on Ban.fiKn.4Fe2As 2 sh ows purely isotropic 
gaps around each surface 57| , which is inconsistent 
with a d-wave model. Although such models have 
been discussed in connection with the iron arsenides 
[2^, most interest has been focussed on extended- 
s± models. In fact, explicit calculations of the neu- 
tron scattering cross section predict the existence 
of a resonant spin excitation with extended-s±, but 
not d-wave symmetry 58 , 5^ . The small Fermi sur- 
faces seen in the iron arsenides would not intersect 
any of the nodal lines in these models, so the gaps 
on each surface would be nearly constant. However, 
the prediction is that the sign of the energy gap 
on the hole pockets at the F-point and the electron 
pockets at the M-point would be opposite, reflecting 
an electron-electron repulsion at short-range but an 
attractive interaction for electrons on neighboring 
iron atoms. The wavevector connecting the hole 
and electron pockets is precisely where the resonant 
spin excitation has been observed. 

More recently, resonant spin excitations have also 
been seen in single crystals of BaFei.84Coo.i6As2 us- 
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ing both pulsed source and triple-axis instruments 
at Oak Ridge [3| ■ This has confirmed that the ex- 
citation is centered at Q=(i,^) within the plane, 
but it is two-dimensional, with intensity spread out 
along Q=Z following a single-ion Fe^+ form-factor. 
Since Tc = 22 K, which is somewhat lower than in 
Bao.6Ko.4Fe2As2, the energy of the resonance is also 
lower at 9.6 meV, in agreement with the previous 
scaling. 

There is a third report of a spin resonance in 
BaFe1.gNio.1A2 (Tc = 20 K) [16], again observed 
at Q=(i,i,Z). Chi et al observe an energy disper- 
sion of about 2meV; the resonance is at 9.1 meV 
at / = but at 7.0 meV at I = ±1. They at- 
tribute this to antiferromagnetic coupling between 
the FeAs layers, although it could also reflect some 
three-dimensional modulation of the Fermi sur- 
faces. Although this work shows evidence of three- 
dimensionality in the superconducting order, it is 
still much more two-dimensional than the spin wave 
excitations of the parent compounds. 

5. Conclusions 

It is far too soon since the discovery of these fas- 
cinating compounds to declare that any of the im- 
portant issues surrounding their superconductivity 
have been settled. As shown in this review, there 
are many unanswered questions concerning the ori- 
gin of the magnetic interactions, the strength of the 
electron-phonon coupling, or the dimensionality of 
the superconducting order parameter. Neverthe- 
less, it is remarkable how much progress has been 
made in such a short time. This reflects the experi- 
ence that the scientific community has gained over 
the past thirty years of studying unconventional su- 
perconductors, first the heavy fermions and then 
the copper oxide superconductors. Many of the in- 
sights gained in those investigations have been di- 
rectly applied to the new iron arsenide supercon- 
ductors, and, in some cases, theories considered but 
then rejected in other systems have become useful 
here. 

This is particularly true for neutron scattering. 
It took several years for the significance of the res- 
onant spin excitations in the copper oxides to be 
appreciated. Now our theoretical understanding 
is well advanced, particularly now that they have 
been seen in other unconventional superconductors. 
This has allowed the new measurements on the iron 
arsenides to be incorporated into theories of the 
superconductivity much more rapidly than before. 



Although the new methods of detecting extended 
s^-wave symmetry have been proposed [6(|, inelas- 
tic neutron scattering is until now the only phase- 
sensitive technique to provide direct evidence in 
these compounds. 

In spite of the many similarities to other su- 
perconductors, the iron arsenides are not carbon 
copies. If, as seems likely, the superconducting pair- 
ing is mediated by spin fiuctuations, they are much 
more itinerant than in the copper oxides or the 
heavy fermions and the superconducting symme- 
try appears to be quite different. The existence of 
the resonance at the same wavevector as the SDW 
does lend some support to the idea that the transi- 
tion from antiferromagnetism to superconductivity 
occurs when the changes in the Fermi surface with 
doping or pressure suppress the nesting instability. 
Although the conditions for magnetic order are no 
longer satisfied, the susceptibility is still sufficiently 
enhanced to favor spin fiuctuation-mediated super- 
conductivity as discussed by Singh in his review 
[3]. The fact that the wavevectors characterizing 
both the spin density wave and superconductivity 
are the same is therefore no coincidence. 

Inelastic neutron scattering experiments are only 
just beginning in these systems. The improvements 
in sample quality and, particularly important for 
neutrons, size will allow more detailed studies over 
a wider range of wavevectors and energies, with sin- 
gle crystal results superceding the earlier polycrys- 
talline data, as has already begun to happen. We 
have also not addressed the many materials science 
issues resulting from phase separation, both intrin- 
sic and extrinsic. Nevertheless, progress is likely 
to be rapid if the experience of the past year is an 
accurate guide. 
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